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Monte Carlo calculations using the codes PENELOPE and GEANT4 have been performed to 
characterize the dosimetric parameters of the new 20 mm long catheter based '^^P beta source 
manufactured by Guidant Corporation. The dose distribution along the transverse axis and the two 
dimensional dose rate table have been calculated. Also, the dose rate at the reference point, the radial 
dose function and the anisotropy function were evaluated according to the adapted TG-60 formalism 
for cylindrical sources. PENELOPE and GEANT4 codes were first verified against previous results 
corresponding to the old 27 mm Guidant ^^P beta source. The dose rate at the reference point 
for the unsheathed 27 mm source in water was calculated to be 0.215 ± 0.001 cGy s~^ mCi~^, for 
PENELOPE, and 0.2312 ±0.0008 cGy s"^ mCi"\ for GEANT4. For the unsheathed 20 mm source 
these values were 0.2908 ± 0.0009 cGy s"^ mCi"^ and 0.311 ± 0.001 cGy s"^ mCi"\ respectively. 
Also, a comparison with the limited data available on this new source is shown. We found non 
negligible differences between the results obtained with PENELOPE and GEANT4. 

I. INTRODUCTION 

One of the major limitations of coronary interventions lies in the appearance of restenosis of the treated vessel. 
Restenosis occurs in 35%-40% of the treated patients in case of angioplasty, whilst the percentage reduces to 22%-32% 
in case coronary stents are used [21 . 

Intravascular brachytherapy appears to be a useful tool to solve such problems 0. PREVENT % and INHIBIT 
01 studies have pointed out the ability of ionizing radiation to inhibit restenosis. Though either gamma and beta 
radiation have been used, endovascular beta-radiotherapy has become most popular in Europe due to the obvious 
radiation safety advantages. 

This new, rapidly growing interest has led both the AAPM Q and the ESTRO to report recommendations on the 
medical physic issues. At the minimum standard level, dose at the prescription point in the central plane is encouraged 
to be recorded and reported. Additional dose calculations in peripheral planes require complete two-dimensional (2D) 
dosimetry information. However, beta sources present a high dose gradient, the distances in which the dose is deposited 
being considerably small. This has made it difficult to perform adequate dosimetry on intravascular sources. 

Recently, the dosimetric characterization of a 27 mm long, catheter based, "^^P beta source wire, manufactured by 
Guidant Corporation (Houston, Texas, USA), has been done by means of plastic scintillator and radiochromic film 

US. 

Monte Carlo (MC) simulations have proved to be a valuable complementary means for determining the dose 
distributions around beta emitters. The ^^P source mentioned above has been simulated with different MC codes 
(MCNPv4B2 0-II3, CYLTRAN/ITS3 0,[I3, MCNPv4C jg, EGS4 and EGSnrc fll) and the calculated 

doses agree within < 10% with the experimental results. 

In this work we address the 2D dosimetric characterization of a new, clinically used, 20 mm long catheter based 
■^^P beta source wire, commercially available in the Galileo"^'^ Intravascular Radiotherapy System, manufactured also 
by Guidant Corporation. This source has been used in a detailed study in which a variety of MC codes are 
intercompared. Also a comparison with measurements done by using radiochromic-dye film and an extrapolation 
chamber is carried out. Measured and calculated values were found to agree within 10% in a polystyrene phantom. 
On the other hand, the calculated depth dose curves, using the MC codes, were within 5% at 4 mm depth, after an 
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FIG. 1: Schematic view of the 27 (upper panel) and 20 mm (lower panel) ^^P sources manufactured by Guidant Corporation. 



adequate selection of the tracking parameters. In this last case water was considered the medium surrounding the 
source. 

Here we have used the MC codes PENELOPE T?] and GEANT4 il6|]. PENELOPE was designed specifically 
to account for the transport of low-energy electrons, photons and positrons and has been used to simulate other 
brachytherapy sources GEANT4 includes special packages which permit to describe low energy physics 

processes and, at present, it is being used to simulate different situations of medical physics interest p^ . In this work 
we have determined the dose rate at the reference point, the depth dose and the two dimensional dose distribution. 
We also calculated the radial dose and the anisotropy functions, as proposed in the adapted AAPM TG-60 formalism 
for cylindrical sources , used also in Ref . ■ We compare some of our results with those quoted in Ref . . 

II. MATERIAL AND METHODS 
A. Source models 

A schematic view of the geometry of the two sources analyzed is shown in Fig. 1. The 27 mm source has been 
described in Refs. 00. The geometry of the 20 mm source has been obtained from Ref. ^md corresponds closely 
to that quoted in The active cores are cylinders with length of 27 and 20 mm, respectively. These cores are 

made of C2He02. The 20 mm source presents a 0.019 mm thick source carrier wrapping the active core. This is made 
of polymide. The cores are encapsulated by NiTi tubes with thickness of 0.076 mm in both sources. In the 27 mm 
source, a tungsten X-ray marker of 1 mm in length follows the core at the distal end. In the case of the 20 mm source, 
two markers of the same characteristics are located at both the distal and the proximal ends of the active core. Both 
sources include a NiTi plug of 1 mm in length following the distal tungsten X-ray marker. A NiTi hemispherical seal 
is welded at the distal end and an unbound NiTi wire is attached to the proximal end in both sources. Empty regions 
between the different pieces and the capsules are supposed to be filled with air. 
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Sources are centered inside a catheter of Nylon-12. In some of the calculations done for the 27 mm source, we have 
considered a helicoidal catheter (labeled catheter 1) reported in Refs. 0-E3- °f 20 mm source, we 

have calculated using both this helicoidal catheter and a lobed catheter (labeled catheter 2) currently in use. This 
catheter is a new design to increase the blood flow during the procedure. Both catheters were modeled as tubes with 
an external radius of 0.4445 mm, while the inner diameter was 0.597 mm for catheter 1 and 0.5716 mm for catheter 
2. Gaps between catheters and sources are assumed to be filled with air. 

Catheter and source were placed in a water phantom the dimensions of which are large enough to ensure no 
significant loss of radiation. 

The composition and density considered in the simulation for the different materials is given in Tabled 





C2H6O2 


NiTi 


Tungsten 


Polymide 


Air 


Nylon-12 


H 


0.097432 






0.026362 




0.11749 


C 


0.387026 






0.691133 


0.000124 


0.73045 


N 








0.073270 


0.755267 


0.07098 


O 


0.515542 






0.209235 


0.231781 


0.08108 


Ar 










0.012827 




Ti 




0.44 










Ni 




0.56 










W 






1.0 








density [g cm~"^] 


1.1155 


6.5 


19.3 


1.096 


0.0012048 


1.196 



TABLE I: Composition of the different materials used in the MC simulations performed in this work. The values correspond 
to the weight fraction of each element in the material. In the case of C2II6O2, the density quoted is that used in the 20 mm 
source. For the 27 mm source a density of 1.12 g cm~^ has been used. 



B. Monte Carlo calculations 

In this work the PENELOPE J^J and GEANT4 [3 Monte Carlo codes have been used to perform the simulations. 
PENELOPE is a general purpose MC code which allows to simulate the coupled electron-photon transport. It can be 
applied to energies ranging from a few hundred eV up to 1 GeV and with arbitrary materials. Besides, PENELOPE 
permits a good description of the particle transport at interfaces and presents an accurate description of the electron 
transport at low energies. These characteristics make PENELOPE a useful tool for medical physics applications as 
previous works have pointed out (see e.g. Refs. [itI fisll') . We have used the version 2001 of the code. Details about 
the physical processes included can be found in Ref. p^ . 

GEANT4 is a toolkit aiming to simulate the passage of particles through matter for applications in high energy 
physics, nuclear experiments, medical physics studies, etc. Due to its purpose, the available energy range and the 
possible particles to be simulated are much larger than for PENELOPE. GEANT4 includes low energy packages which 
extend the energy range of particle interactions for electrons, positrons and photons down to 250 eV, and can be used 
up to approximately 100 GeV. In this work we have used the version 4.1 with the patch-01 and the low energy package 
G4EML0W1.1. Detailed documentation concerning the physics involved can be found in Ref. [T6| . 

As the sources considered here present cylindrical symmetry, a cylindrical coordinate system was used with the 
origin in the center of the active core and with the z axis along the source and pointing to its distal end. The dose 
distributions were supposed to be functions of the axial coordinate z and the radial coordinate p. The position of the 
reference point is (po — 2 mm, zq = 0). 

The radionuclide was assumed to be uniformly distributed inside the active core. The direction of the emitted 
electrons was supposed to be isotropic around the initial position. The initial electron energies were sampled from the 
corresponding Fermi distributions, taken from Ref. 2 0J. A calculated average energy of 0.6957 MeV was obtained, 
in agreement with values reported in previous works [a, l2l| . 

In PENELOPE analog simulation is performed for photons. For electrons, instead, the simulation is done in a mixed 
scheme in which collisions characterized by a scattering angle larger than a certain value are called hard collisions 
and are individually simulated, while collisions with a scattering angle smaller that the limit value are called soft 
collisions and are described by means of a multiple scattering theory. The electron tracking is controlled by means of 
four parameters. Ci and C2 refer to elastic collisions. Ci gives the average angular deflection due to a elastic hard 
collision and to the soft collisions previous to it. C2 represents the maximum value permitted for the average fractional 
energy loss in a step. On the other hand, Wcc and Wa are energy cutoffs to distinguish hard and soft events. Thus, 
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the inelastic electron collisions with energy loss W < Wcc and the emission of Bremsstrahlung photons with energy 
W < Wcr are considered in the simulation as soft interactions. In our calculations the simulation parameters were 
fixed to the following values: Wcc = 5 keV, Wcr = 1 keV, Ci = C2 = 0.05, Smax = lO'^^. This last parameter is an 
upper bound for the length of the steps generated in the simulation. 

Photons were simulated down to 1 keV. Primary electrons, S rays and positrons were absorbed when they slow 
down to kinetic energies of 1 keV, in air, and 50 keV, in all the remaining materials, except in the active core where 
the cut off energy was fixed to 200 keV. This value takes care of the fact that electrons with such an energy are 
not expected to leave the active core. Those betas generated with energies below the simulation threshold were not 
transported and their energy was deposited at the source position. 

In case of GEANT4 simulations, it is necessary to choose a given range threshold for each particle type. This 
threshold is converted, internally, to an energy threshold below which secondary particles are not emitted. In our 
calculations the range thresholds have been fixed to 100 nm providing energy thresholds of 1 kcV for both photons 
and electrons in all materials considered. 

The statistical uncertainties were calculated by scoring both the energy deposited in each voxel and its square for 
each history. These voxels were taken to be annular bins with thicknesses Ap ~ 0.1 mm and Az = 0.2 mm. The 
average energy deposited in a given voxel (per incident particle) is 

1 ^ 
E - — Ve 

i=l 

where N is the number of simulated histories and is the energy deposited by all the particles of the i-th history 
(that is, including the primary particle and all the secondaries it generates). The statistical uncertainty is given by 



\ 



1 

N 



1 ^ 



The uncertainties given throughout the paper correspond to Icr. In tables, and in order to simplify the writing, they 
are given between parentheses; for example, a value quoted as 0.956(2) means 0.956 ± 0.002. 

In case of the 27 mm source, the simulations have been done by following 2 x 10^ histories for PENELOPE and 
5 X 10^ for GEANT4. For the 20 mm source a total of 5 x 10^ histories were simulated with both codes. These 
numbers permitted to keep the statistical uncertainties under reasonable levels. 

Additionally, we performed, for the 20 mm source, simulations with PENELOPE by using scoring voxels given by 
cylindrical annuli of Az — 13 mm centered along the source axial direction, and Ap = 0.1 mm in the radial direction. 
A total of 6 x 10^ histories were followed in this case. These scoring voxels warranty radial distributions which are 
effectively constant over the axial length and permit the uncertainty to be reduced considerably. 



III. RESULTS 



A. 27 mm source 



First we validated the MC codes with the old 27 mm source. Our results were compared with those reported in 
literature. 

In Refs. and the MCNPv4B2 code was used. Mourtada et al. transported 10^ histories in each run, with 
scoring voxels given by cylindrical annuli of 20 mm in the source axial direction and 0.1 mm in the radial direction. 
They showed also the dose rate per unit activity along the transverse axis. The same code and voxels were used by 
Bohm et al. |1] for the dose rate per unit activity along the transverse axis. Besides, these authors calculated the dose 
to water, away and along the source, using thin cylindrical shells with a thickness of 0.1 mm and a height of 0.1 mm 
as scoring voxels. Wang and Li 12] used scoring voxels equal to those we have considered here. Finally, Seltzer 
quoted the results obtained for a simulation with CYLTRAN/ITS3, but did not give information about the scoring 
voxels. 

The dose rate at the reference point for the unsheathed 27 mm source in water was calculated to be 0.215 ± 0.001 
and 0.2312 ± 0.0008 cGy s^^ mCr^ for PENELOPE and GEANT4, respectively, thus differing by 7%. These 
values must be compared with the results obtained with the other codes. The calculation with EGSnrc 1^ provides 
0.2185 ±0.0002 cGy s'^ mCi'^ Mourtada et al. obtained 0.229 cGy s~i mCi"! by means of MCNPv4B2. Finally, 
in Ref. a value of 0.232 ±0.00028 cGy s'^ mCi"! was found also with MCNPv4B2. This is the largest value found, 
which differs from the one we have obtained with PENELOPE (the smallest one) by ~ 8%. 
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z 


P 


this work 


Ref. [8J 


Ref. fll] Ref. 


fl2] 


[mm] 


[mm] 


PENELOPE GEANT4 MCNPv4B2 CYLTRAN/ITS3 EGSnrc 


EGS4 


0.0 


1.0 


0.718(3) 


0.721(2) 


0.728 


0.801 0.731 


0.729 




2.0 


0.221(1) 


0.2379(9) 


0.232 


0.232 0.225 


0.222 




3.0 


0.0718(6) 


0.0868(4) 


0.0785 


0.0767 0.0732 


0.0708 




4.0 


0.0202(3) 


0.0281(2) 


0.0229 


0.0212 0.0207 


0.0193 


10.0 


1.0 


0.719(3) 


0.717(2) 


0.725 








2.0 


0.220(1) 


0.2356(9) 


0.230 








3.0 


0.0708(6) 


0.0853(4) 


0.0769 








4.0 


0.0206(3) 


0.0275(2) 


0.0223 






12.0 


1.0 


0.677(3) 


0.667(2) 


0.681 








2.0 


0.193(1) 


0.2046(8) 


0.202 








3.0 


0.0609(6) 


0.0713(4) 


0.0654 








4.0 


0.0167(3) 


0.0230(2) 


0.0186 






14.0 


1.0 


0.180(2) 


0.187(1) 


0.188 








2.0 


0.0755(8) 


0.0816(5) 


0.0793 








3.0 


0.0266(4) 


0.0322(3) 


0.0291 








4.0 


0.0076(2) 


0.0109(1) 


0.00870 






16.0 


1.0 


0.0097(5) 


0.0130(3) 


0.0113 








2.0 


0.0067(3) 


0.0101(2) 


0.00858 








3.0 


0.0035(1) 


0.0051(1) 


0.00423 








4.0 


0.00111(7) 


0.00192(6) 


0.00135 







TABLE II: Two-dimensional dose rates (in cGy s~^ mCi~^) for the 27 mm source sheathed by catheter 1. The errors corre- 
sponding to the calculations performed with MCNPv4B2 (Ref. j^j) as well as those with EGSnrc and EGS4 (Ref. ^'^^ 
quoted to be below 1%. The CYLTRAN/ITS3 data have been obtained by scanning Fig. 11 of Ref. [Tll | and therefore are 
indicative only. 

In Table |n] we compare the dose rate per unit activity along the transverse axis (values for z = in the first set 
of rows) obtained with PENELOPE and GEANT4 with those of Refs. |l|ll|l3- The resuhs correspond to the 
wire source sheathed by catheter 1. As we can see PENELOPE produces results very similar to those of EGSnrc. 
In this case the differences for any p value are below 6% with respect to the dose rate obtained with PENELOPE 
at the reference point, -Dpen(poj ^o)- Up to p = 3 mm PENELOPE practically coincides with EGS4, the differences 
relative to -Dpen(/Oo, ^o) being smaller than 5%. The difference between GEANT4 (MCNPv4B2) and PENELOPE 
is below 8% (5%), relative to I)pen(poj ^o)- In what respect to the results of CYLTRAN/ITS3, they show negligible 
differences with PENELOPE results. The fact that these last quoted values have been obtained by scanning Figure 
11 of Ref. makes these differences to be only indicative. 

Two-dimensional dose distributions are shown in Tables HTl and UTTl for some values of p and z > 0. Similar results 
are obtained for the negative values of z, though the source wire is not symmetric with respect to z = 0. In Table Hll 
our results are compared with those of Bohm et al. calculated with MCNPv4B2 for the source sheathed by catheter 
1. In Table Unl the dose distributions for the unsheathed source, normalized to the reference point, are compared with 
the results obtained by Wang and Li jl^] using EGSnrc. 

The differences between MCNPv4B2 and PENELOPE (see TableHHl, remain below 5%, relative to I)pen(po, ^q)). 
Same situation appears for MCNPv4B2 and GEANT4, which give very similar results for p < 2 mm. 

Table ITlll shows that EGSnrc and PENELOPE provide very similar normalized doses. It is important to remember 
here that the difference in the doses at the reference point between both codes is 1.6% only. The differences maintain 
at more or less the same level for all z and p values studied. 

Finally, and as we can see in both tables, the discrepancies between PENELOPE and GEANT4, observed for the 
dose rate per unit activity along the transverse axis, show up again. In any case the differences remain below 8% with 
respect to Dpen(po,zo). 
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p 








z [mm] 






[mm] 


MC code 


0.0 


10.0 


12.0 


14.0 


16.0 




EGSmc 


3.2844 


3.2730 


3.0795 


0.8230 


0.0439 


1.0 


PENELOPE 


3.31(1) 


3.29(1) 


3.11(1) 


0.835(8) 


0.062(2) 




GEANT4 


3.036(9) 


3.023(9) 


2.835(8) 


0.789(5) 


0.051(1) 




EGSmc 


1.0000 


0.9926 


0.8762 


0.3363 


0.0324 


2.0 


PENELOPE 


1.000(6) 


0.996(6) 


0.881(6) 


0.336(3) 


0.039(1) 




GEANT4 


1.000(4) 


0.988(4) 


0.860(3) 


0.346(2) 


0.0394(8) 




EGSmc 


0.3226 


0.3194 


0.2711 


0.1176 


0.0158 


3.0 


PENELOPE 


0.323(3) 


0.317(3) 


0.271(3) 


0.117(2) 


0.0182(7) 




GEANT4 


0.358(2) 


0.352(2) 


0.299(2) 


0.131(1) 


0.0209(5) 




EGSmc 


0.0897 


0.0890 


0.0744 


0.0336 


0.0050 


4.0 


PENELOPE 


0.089(1) 


0.091(1) 


0.074(1) 


0.0329(8) 0.0050(3) 




GEANT4 


0.1153(9) 0.1125(9) 0.0957(8) 0.0453(6) 0.0079(2) 



TABLE III: Two-dimensional dose rates for the unsheathed 27 mm source. The values are normalized to the dose rate at the 
reference point (po ~ 2 mm, zo = 0). The errors corresponding to the calculations performed with EGSnrc (Ref. are 
quoted to be below 1%. 



B. 20 mm source 



1. Unsheathed source 



Once we have analyzed the resuhs obtained for the 27 mm source, we wih discuss the dose rates found for the 
20 mm source. AU the results here correspond to the unsheathed source. For clarity, we leave the analysis of the 
sheathed source for the next subsection. 

In Table IIVI (box labeled z = 0) we show the results we have obtained for the dose rate per unit activity along 
the transverse axis. The first row gives the results obtained using PENELOPE for voxels with Az — 13 mm. The 
second (third) rows correspond to the results obtained with PENELOPE (GEANT4) for voxels with Az = 0.2 mm. 
The dose rates at the reference point calculated with PENELOPE are 0.2918 ± 0.0001 cGy s^^ mCi~^ when the large 
scoring voxel (that with Az — 13 mm) is used and 0.2908 ± 0.0009 cGy s~^ mCi~^ for the small one (that with 
Az = 0.2 mm ). GEANT4 produces a value of 0.311 ± 0.001 cGy s'^ mCi-\ which differs from the previous ones m 
~ 7%. On the other hand, the discrepancy between the results obtained with both codes increases with increasing p, 
showing a behavior similar to that found for the 27 mm source. 




FIG. 2: Axial dose distributions for the unsheathed 20 mm source in water, for p = 1 to 6 mm. The values are normalized to 
the reference point (po = 2 mm, zq = 0). PENELOPE results are plotted with points, while curves correspond to GEANT4. 
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FIG. 3: Radial dose distributions for the unsheathed 20 mm source in water, for 2 = and 2; = 7 to 14 mm. The values 
are normalized to the reference point (po = 2 mm, zo = 0). PENELOPE results axe plotted with points. GEANT4 results 
corresponding to 2: = 0, 11 and 14 mm are plotted with lines. 
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FIG. 4: Ratio of the dose at different z values to the dose at —z as a function of p, for the unsheathed 20 mm source in water. 
PENELOPE (black squares) and GEANT4 (open circles) results are shown. 



at — z. The results are similar in both codes. As we can see, the ratio of the dose rates does not show statistically 
significant differences from unity, even at large z. This means that the asymmetry of the source structure due to the 
small differences between the distal and the proximal ends is smaller than statistical uncertainties and would therefore 
probably be irrelevant in practice. 

In Fig. 5 we analyze the axial dose uniformity by plotting the ratio of the dose at z = 3 up to 8 mm to the dose 
at 2; = as a function of p. The results show that, for a given p value, doses stay constant for axial positions up 
to 6-7 mm. This occurs for both PENELOPE and GEANT4 and means that the scoring voxels with Az — 13 mm, 
considered above to reduce the uncertainty in calculating the radial dose, can be used safely for that purpose. 

The two-dimensional dose rates we have obtained using voxels with Az = 0.2 mm are listed in Table Hvl Eor each 
z value, the first row corresponds to PENELOPE and the second one to GEANT4. The values for z < are not 
included because of the symmetry with respect to the corresponding z > results. On the other hand, the uniformity 
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+ PENELOPE GEANT4 




p [mm] p [mm] 



FIG. 5: Ratio of the dose at different z values to the dose at 2 = as a function of p, for the unsheathed 20 mm source in 
water. Results obtained with PENELOPE (GEANT4) are shown with black squares (open circles). 



up to z 6 mm is again evident. The comparison between both codes indicates that GEANT4 predicts larger doses 
than PENELOPE except for p below 0.5 mm. 

To finish with the analysis of the unsheathed source, we calculated the radial dose functions and the anisotropy 
functions. Following Ref. |l9l |. these quantities are given in the cylindrical coordinate system instead of the spherical 
one used in the TG-60 formalism ^ . For the radial dose we have 

D{p,zq) G{po,zo) 
9{P) = r 7 ■ 



p 

mm] 


20 mm 


source 


27 mm 


source 


PENELOPE GEANT4 


PENELOPE EGSnrc 


0.5 


1.698(7) 


1.507(6) 


1.75(1) 


1.7410 


1.0 


1.542(7) 


1.428(6) 


1.57(1) 


1.5582 


1.5 


1.275(6) 


1.226(6) 


1.30(1) 


1.2867 


2.0 


1.000(5) 


1.000(5) 


1.000(9) 


1.0000 


2.5 


0.742(4) 


0.780(4) 


0.730(7) 


0.7343 


3.0 


0.518(3) 


0.586(3) 


0.511(5) 


0.5093 


3.5 


0.342(2) 


0.412(3) 


0.334(4) 


0.3310 


4.0 


0.208(2) 


0.270(2) 


0.199(3) 


0.1990 


4.5 


0.119(1) 


0.167(2) 


0.108(2) 


0.1098 


5.0 


0.059(1) 


0.090 (1) 


0.054(2) 


0.0543 


5.5 


0.0257(6) 


0.0451(9) 


0.025(1) 


0.0231 


6.0 


0.0099(4) 


0.0190(6) 


0.0093(7) 


0.0083 


6.5 


0.0028(2) 


0.0070(3) 


0.0034(4) 


0.0025 


7.0 


0.0006(1) 


0.0017(2) 


0.0008(2) 


0.0007 



TABLE V: Radial dose functions for the P sources. Simulations were performed in water. 
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FIG. 6: Comparison between the anisotropy functions for the unsheathed 20 mm source obtained with PENELOPE (squares) 
and GEANT4 (curves). The errors for GEANT4 results are of the same order than those obtained with PENELOPE. 



The anisotropy function is given by 

^ D{p,z) G{p,Zo) 

where G{p, z) — l3/{Lp), with L = 27 or 20 mm, according to the source considered, and (3 is the angle between the 
hues joining the point and the two extremes of L. 

Calculated radial dose functions for PENELOPE and GEANT4 in water are shown in Tabled As we can see, the 
discrepancies observed between both codes in previous results are evident here again. Besides, we have included the 
results obtained with PENELOPE for the 27 mm source and those quoted by Wang and Li for EGSnrc (T^. The 
radial dose function for the 20 mm source appears to be rather similar to that of the 27 mm source. The results 
obtained with PENELOPE differs less than 10% up to p = 6 mm. Finally, PENELOPE and EGSnrc provide resuhs 
which coincide within the statistical uncertainties for the 27 mm source. 

Table IVII shows the anisotropy function calculated for the unsheathed 20 mm source with PENELOPE. These 
results are compared in Fig. 6 with those obtained with GEANT4. As in previous comparisons, the discrepancies 
between both codes appear to be noticeable for large p and/or z values, though the absolute differences remain small 
in comparison with the values obtained at the reference point. 



2. Sheathed source 



Here we analyze the effect of the catheter on the dosimetry. First we compare the results we have obtained with both 
codes with the scarce data available in the literature. Table IVTll shows this comparison with the data provided by the 
manufacturer in the source manual (second column) and with the results of a set of calculations with different codes 
performed by Seltzer pj|. In these calculations the scoring voxels were taken to be annuli of radial thickness 0.1 mm 
and axial length 15 mm, centered along the active wire length. This produces greatly reduced statistical uncertainties. 
As can be seen, our results for PENELOPE coincide with those found with CYLTRAN/ITS3 and PENELOPE in 
Ref. except for small deviations (less than 1%) at p = 2 mm in the last case. The values provided by Guidant 
are between 5 and 10% above those we have found for PENELOPE. MCNP4G results are similar to the Guidant ones 
and those calculated with EGS4 and EGSnrc fall to ~ 7% below those obtained with PENELOPE at p = 4 mm. The 
values found with GEANT4 follow the same trend seen in the previous sections. 

Fig. 7 shows the relative difference 
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PENELOPE ^ GEANT4 




1 2 3 4 5 

p [mm] 



FIG. 7: Relative difference between the depth dose rates calculated for the 20 mm source sheathed, Z)cath(p, zo), and unsheathed, 
D{p,zo), with PENELOPE (squares) and GEANT4 (disks). Upper (lower) panel corresponds to catheter 1 (2). 



between the depth doses calculated with and without catheter for PENELOPE (squares) and GEANT4 (disks). The 
results corresponding to the 20 mm source are shown. As we can see, both codes produce similar results with both 
catheters. However, GEANT4 shows an increase of the effect of the catheter above p = 3.5 mm which is not present 

in the PENELOPE calculation. Similar results arc obtained for the 27 mm source sheathed by catheter 1. On the 
other hand, catheter 2 produces a modification of the dose in water slightly smaller than catheter 1. This last gives 
rise to enhancements of the dose rate on the transverse axis reaching ~ 4% above p = 3 mm, whilst the first one 
increases the dose rate by ^ 2.5% at most. This result is reasonable if one considers that catheter 2 has a smaller 
inner radius than catheter 1, thus reducing the air gap. In any case, uncertainties in source strength measurements 



z 

[mm] 










P 


[mm] 










0.4 


0.6 


0.8 


1.0 


1.5 


2.0 


2.5 


3.0 


4.0 


5.0 


0.0 


1.000(3) 


1.000(3) 


1.000(3) 


1.000(4) 


1.000(4) 


1.000(5) 


1.000(6) 


1.000(7) 


1.00(1) 


1.00(2) 


2.0 


1.001(3) 


1.003(3) 


0.999(3) 


0.998(4) 


1.014(4) 


1.007(5) 


1.009(6) 


1.012(7) 


1.00(1) 


1.03(2) 


4.0 


1.005(3) 


1.009(3) 


1.006(3) 


1.009(4) 


1.024(4) 


1.022(5) 


1.035(6) 


1.038(7) 


1.05(1) 


1.02(2) 


6.0 


1.013(3) 


1.022(3) 


1.025(4) 


1.028(4) 


1.055(4) 


1.071(5) 


1.094(6) 


1.091(8) 


1.12(1) 


1.14(3) 


7.0 


1.022(3) 


1.039(3) 


1.042(4) 


1.050(4) 


1.087(4) 


1.111(5) 


1.130(6) 


1.153(8) 


1.16(1) 


1.17(3) 


8.0 


1.038(3) 


1.057(3) 


1.067(4) 


1.084(4) 


1.122(5) 


1.148(5) 


1.151(7) 


1.169(8) 


1.16(1) 


1.14(3) 


9.0 


1.068(3) 


1.094(4) 


1.108(4) 


1.115(4) 


1.128(5) 


1.128(6) 


1.102(7) 


1.106(8) 


1.08(1) 


1.05(3) 


10.0 


0.990(4) 


0.977(4) 


0.967(4) 


0.949(4) 


0.944(5) 


0.927(5) 


0.905(6) 


0.894(7) 


0.87(1) 


0.81(2) 


11.0 


0.261(4) 


0.331(4) 


0.392(4) 


0.441(4) 


0.512(4) 


0.548(5) 


0.561(5) 


0.569(6) 


0.542(9) 


0.49(2) 


12.0 


0.128(4) 


0.153(3) 


0.168(3) 


0.179(3) 


0.218(3) 


0.246(4) 


0.260(4) 


0.269(4) 


0.250(6) 


0.24(1) 


13.0 


0.056(3) 


0.057(2) 


0.061(2) 


0.067(2) 


0.077(2) 


0.086(2) 


0.089(2) 


0.095(3) 


0.084(4) 


0.079(5) 


14.0 


0.016(2) 


0.018(2) 


0.018(1) 


0.018(1) 


0.018(1) 


0.023(1) 


0.023(1) 


0.026(2) 


0.019(3) 


0.017(6) 


15.0 


0.005(1) 


0.0025(7) 0.0037(9) 0.0018(4) 0.0034(7) 0.0038(6) 0.0031(4) 0.0032(6) 0.00351(3) 0.0068(1) 



TABLE VI: Anisotropy function F{p, z) for the unsheathed 20 mm source in water, calculated with PENELOPE. 
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p 

mm] 




Ref. [1^ 


this work 


Guidant CYLTRAN/ITS3 MCNP4C EGS4 EGSnrc PENELOPE 


PENELOPE GEANT4 


1.0 


0.999 




0.975(2) 


0.971 (2) 


1.5 


0.548 




0.526(1) 


0.540 (1) 


2.0 


0.317 


0.300 0.314 0.296 0.296 0.298 


0.3005(9) 


0.321 (1) 


2.5 


0.184 




0.1723(6) 


0.1957(7) 


3.0 


0.105 




0.0976(4) 


0.1172(5) 


3.5 


0.0569 




0.0534(3) 


0.0683(4) 


4.0 


0.0290 


0.0276 0.0285 0.0255 0.0254 0.0274 


0.0274(2) 


0.0384(3) 


4.5 


0.0135 




0.0133(1) 


0.0200(2) 



TABLE VIL Depth doses (in cGy s ^ mCi ^) for the 20 mm source sheathed by catheter 1. The values of the second column 
are those provided by the manufacturer in the GALILEO"""^ system manuaL Simulations were performed in water. 



might outweigh this correction in chnical practice. 



IV. CONCLUSIONS 



In this work we have performed the dosimetric characterization of two wire sources of '^^P (of 27 mm and 20 mm 
active length) manufactured by Guidant Corporation and clinically used in intravascular brachytherapy. To do so we 
have used the Monte Carlo codes PENELOPE and GEANT4. 

First we have validated the codes by comparing the results obtained with these two codes with those available in 
the literature for the 27 mm source. The main conclusion is that PENELOPE produces dose rates in reasonable 
agreement with those quoted in Ref. for EGSnrc and EGS4 except in the tails of the axial doses (large z values). 
Calculations performed with MCNPv4B2 0, |1| differ from those obtained with PENELOPE as the distance to the 
source increases. The comparison between PENELOPE and GEANT4 shows a clear disagreement. 

For the 20 mm source, we have calculated the depth doses on the transverse axis, the two-dimensional dose rate 
table and the radial dose and anisotropy functions (this two last according to the approach proposed in Ref. 
for cylindrical sources). We have checked that the dose distribution remains uniform along the axial direction for 
\z\ smaller than 6-7 mm and p 7 mm. Besides, we have tested the symmetry of the dose distribution. The dose 
rates calculated with PENELOPE show a good agreement with the scarce data available, while those obtained with 
GEANT4 follow a trend similar to that observed in the 27 mm source. 

The presence of the catheter sheathing the sources produces an increase in the dose rate on the transverse axis. 
This enhancement is larger for the helicoidal catheter (catheter 1) than for the lobed one (catheter 2). For the first 
one, the dose on the transverse axis increases up to 4% above p — 3 mm, while the second one produces a maximum 
enhancement of ~ 2.5%. These differences could produce non-negligible modifications in the clinical use of the sources. 

The results we have obtained indicate that PENELOPE produces a dose rate map which can be considered to rep- 
resent the dosimetric characterization of the sources here studied. The appreciable discrepancies shown by GEANT4 
results need a deeper investigation in order to clarify the origin of these differences. In particular, to determine if 
they are due to the physics included in the simulations or to the tracking procedure should be of great interest. We 
are working in this direction. 

The full data set we have obtained here is available at http://fml31.ugr.es/P32-IVB 
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